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Abstract 

It is studied the possibility of direct extraction of the transversity and its accom¬ 
panying T-odd parton distribution function (PDF) from Drell-Yan (DY) processes with 
unpolarized pion beam and with both unpolarized and transversely polarized proton tar¬ 
gets. At present, such a measurement can be performed on the COMPASS experiment at 
CERN. The preliminary estimations performed for COMPASS kinematic region demon¬ 
strate that it is quite real to extract both transversity and its accompanying T-odd PDF 
in the COMPASS conditions. 


The advantage of DY process for extraction of PDFs, is that there is no need in any frag¬ 
mentation functions. It is well known that the double transversely polarized DY process 
H\h\ —> l + l~X allows to directly extract the transversity distributions (see Ref. [Ij for review). 
In particular, the double polarized DY process with antiproton beam is planned to be studied at 
PAX |2j. However, this is a rather difficult task to produce antiproton beam with the sufficiently 
high degree of polarization. So, it is certainly desirable to have an alternative (complementary) 
possibility allowing to extract the transversity PDF from unpolarized and single-polarized DY 
processes. This could be a matter of especial interest for the COMPASS experiment p j where 
the possibility to study DY processes with unpolarized pion beam and with both unpolarized 
and transversely polarized proton targets ir~p —*• p + p~X, tt~ p' 1 —> p + p~X is under discussion 
now. 

The original expressions for unpolarized and single-polarized DY cross-sections |3j are very 
inconvenient in application since all An-dependent PDFs enter there in the complex convolution. 
To avoid this problem in Ref. |5] the qr integration approach P0IB! was applied. As a result, 
the procedure proposed in Ref. [5] allows to extract the transversity hi and the first moment 


hif\x) = / d 2 k 


T {ml) Afe.R) 


( 1 ) 


of T-odd distribution hf directly, without any model assumptions about fc^-dependence of 
hf(x, kji). 

The general procedure proposed in Ref. |5] applied to unpolarized DY process tt~p —> 
p + p~X gives 1 
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'Eq. m is obtained within the quark parton model. It is of importance that the large values of v cannot 
be explained by leading and next-to-leading order perturbative QCD corrections as well as by the high twists 
effects (see @| and references therein). 


1 










where k is the coefficient at cos 20 dependent part of the properly integrated over qt ratio of 
unpolarized cross-sections: 
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At the same time in the case of single polarized DY process tt pi —> /i + p X, operating just as 
in Ref. [£], one gets 


1 T, q e 2 q [h^ 1 \x 7r )h lq (x p ) + (x n <-»• x p )] 
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2 T, q e 2 q [fi q (x n )f lq (x p ) + (x n <-+ x p )} 
where the single spin asymmetry (SSA) Ah is defined as 2 
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In Eqs. (EH) the quantity h 


i(i), 

1 q 
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r ) is defined by Eq. ©• All other notations are the same as 


in Ref. 0 (see Ref. [T] for details on kinematics in the Collins-Soper frame we deal with). 
Neglecting strange quark PDF contributions, squared sea contributions of w-quark PDF 
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, fiu(%n) _fiu{x P ) , and cross terms containing the products of sea and 


valence d-quark PDFs (additionally suppressed by the charge factor 1/4), one arrives at the 
simplified equations 
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Notice that while two equations corresponding to unpolarized and single-polarized antiproton 
proton DY processes completely determine the transversity and its accompanying T-odd PD 


in proton 0, two Eqs. m and ( 0 ) contain three unknown quantities h ,/ 7 x 7 ) _, h u ) (x p ) 

7 r 

and hi u (xp). Nevertheless, from Eqs. CD and (0) it immediately follows that 


h 




x p ) 


h\ u (xp) 


1 k(xn,x p ) 

7r p 

1 6A h (x„,x p ) 

7 T~p^ 


(9) 


Thus, using only unpolarized pion beam colliding with unpolarized and transversely polar¬ 
ized protons, it is possible to extract the ratio / h lu (x p ) . However, it is certainly 


desirable to extract h^\x p ) 


and hi u (xp) in separation. 

P P 


~ sin(0— ef>s) -ft p — 

^Notice that SSA Ah is analogous to asymmetry A UT N (weighted with sin(0 — 0g) and the same 


weight qr/Misr) applied in Ref. jj| with respect to Sivers function extraction from the single-polarized DY 
processes. 
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The simplest way to solve this problem is to use in Eqs. (0), ® the quantity 

p 

extracted from k measured in unpolarized DY process, pp —> l + l~X (in the way proposed in 
Ref. 0). However, if one wishes to extract all quantities within the experiment with the 
pion beam (COMPASS here) the additional assumptions connecting pion and proton PDFs are 
necessary. Taking into account the probability interpretations of hf q and f\ q PDFs, it is natural 
to write the relation 
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Notice that the assumption given by Eq. m is in accordance (but is much less strong restric¬ 
tion) with the Boer’s model (see Eq. (50) in Ref. [3j), where C u = M p c%/M n c™. 

As we will see below, one should put C u to be about unity 


C u ~ 1 


to reconcile the results on 11 in proton obtained from the simulated k 


7T p 
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with the respective 


results 0 obtained from the simulated k _ as well as with the upper bound |5j on this quantity. 


pp 


Thus, Eqs. © and (JHJ) are rewritten as (c.f. Eqs. (19), (20) in Ref. |5J) 
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Looking at Eqs. dn, m , one can see that now the number of equations is equal to the 
number of variables to be found. Measuring the quantity k in unpolarized DY (Eqs. Q. (©) 
and using Eq. (TH21) one can obtain the quantity h^. V} . Then, measuring SSA, Eq. ©, and 

p 

using in Eq. m the obtained from unpolarized DY quantity , one can eventually extract 

p 

the transversity distribution hi u . 

p 

To deal with Eqs. ca and m in practice, one should consider them at the points 3 
Xn = x p = x (i.e., xf = x n — x p = 0), so that 


and 






7T p 


hi u (x) 





where now all PDFs refer to proton. 

3 The different points xf = 0 can be reached changing Q 2 value at fixed s = Q 2 /x 1 X 2 = Q 2 /r 
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Figure 1: k versus x n for xp ~ 0. Data is obtained with MC simulations for 60 GeV (closed 
circles) and 100 GeV (open circles) pion beams. 


To estimate the possibility of 1 ^ measurement, the special simulation of unpolarized DY 
events with the COMPASS kinematics are performed. The pion-proton collisions are simulated 
with the PYTHIA generator nm- Two samples are prepared corresponding to 60 GeV and 100 
GeV pion beams. Each sample contains about 100 K pure Drcll-Yan events. The events are 
weighted (see Ref. |5J for detail) with the ratio of DY cross-sections given by (see Refs, pll II Iji 


daW/dn 

R = ^ ’ 

3 

R = -(1 + cos 2 6 + (v/2) cos2</>sin 2 0) (v = 2k), 
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where v dependencies of qp and x n are taken from Refs. mm 
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Figure 2: versus x n for xp — 0. Data is obtained with MC simulations for 60 GeV (closed 

circles) and 100 GeV (open circles) pion beams. 


The angular distributions of R (Eqs. Q and ©) for both samples are studied just as it 
was done in Ref. [TT] with respect to R (Eqs. (HTJ1) . (fT71) i. The results are shown in Fig. |T] The 
value of k at averaged Q 2 for both energies are found to be 0.7 ± 0.1 for 60 GeV and 0.9 ± 0.1 
for 100 GeV pion beams. 

The quantity hf^ is reconstructed from the obtained values of k using Eq. (fT-ll) with 
Xp = 0 ± 0.04. The results are shown in Fig. [21 Let us recall that to obtain from k , 


7r p 
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Figure 3: SSA given by Eq. m versus xp for 100 GeV pion beam (Qf a , erage = 6.2 GeV 2 ). 


we chosen C u ~ 1. Notice that namely this choice of C u is consistent with the results on 

(compare Fig. [21 with Fig. 4 of Ref. j^j), and also 


obtained in Ref. 


from simulated k 


with the upper bound on estimated in that paper. Otherwise, if C u essentially differs 

from unity, one should multiply the results on by the factor 1/y/Cf, that would lead to 

disagreement of the results on obtained from the simulated quantities k and k _ . 

7 T~p pp 

Certainly, all conclusions made on the basis of simulations are very preliminary. The reliable 
conclusion about C u can be made only from the future measurements of k for both DY processes 
with p and tt~ participation. 

Using the obtained magnitudes of we estimate the expected SSA given by Eq. 03. 
The results are shown in Figs. Q] and [H For estimation of h\ u entering SSA together with h^ tI l 
(see Eq. G3D) we follow the procedure of Ref. m and use (rather crude) “evolution model” 
PE] , where there is no any estimations of uncertainties. That is why in (purely qualitative) 
figures 01 and 0] we present the solid curves instead of points with error bars. To obtain these 
curves we reproduce x-dependence of in the considered region, using the Boer’s model 
(Eq. (50) in Ref. Q), properly numerically corrected in accordance with the simulation results. 



Figure 4: SSA given by Eq. m versus xp for 60 GeV pion beam (Qa Verage = 5.5 GeV 2 ). 

It should be noticed that the estimations of k and A h magnitudes obtained it this paper 
are very preliminary and show just the order of values of these quantities. For more precise 
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estimations one needs the Monte-Carlo generator more suitable for DY processes studies (see, 
for example. Ref. [75] ) than PYTHIA generator which we used (with the proper weighting of 
events) here. 

In summary, it is shown that the proposed in Ref. 5| procedure can be applied to DY 
processes: tt~ p —> p + p~X and -K~p^ —> p + p~X, which could be studied in the COMPASS 
experiment at CERN. The preliminary estimations for COMPASS kinematical region show the 
possibility to measure both k and SSA Ah and then to extract the quantities and h i we 
are interesting in. 
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